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Abstract: In this paper we study the description of the Universe based on the 
low energy superstring theory modified by the Loop Quantum Gravity effects. 
This approach was proposed by De Risi et al. in the Phys. Rev. D 76 (2007) 
103531. We show that in the contrast with the string motivated pre-Big Bang 
scenario, the cosmological realisation of the i-duality transformation is not nec- 
essary to avoid an initial singularity. In the model considered the universe starts 
its evolution in the vacuum phase at time t — > — oo. In this phase the scale 
factor a — * 0, energy density p — > and coupling of the interactions g\ — > 0. 
After this stage the universe evolves to the non-singular hot Big Bang phase 
P Pmax < oo. Then the standard classical universe emerges. During the whole 
evolution the scale factor increases monotonically. We solve this model analyti- 
cally. We also propose and solve numerically the model with an additional dilaton 
potential in which the universe starts the evolution from the asymptotically free 
vacuum phase g 2 s — > and then evolves non-singularly to the emerging dark 
energy dominated phase with the saturated coupling constant g 2 — > const. 
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1. Introduction 

It is commonly believed and physically motivated that gravity exhibits quantum 
nature close to the Planck energy scale. However it can be realised that such 
an energy scale will be never reached with the present generation of accelerators. 
The fact that there is lack of experimental verification, is the main handicap in 
development of the quantum version of General Relativity. However according 
to the standard Big Bang scenario, our Universe arose from the epoch in which 
quantum gravitational description is adequate. This potentially allows for the 
verification of the quantum gravitational models by the astronomical observa- 
tions. On the other side the application of the quantum gravity effects in the 
description of the universe could help to answer the most fundamental questions 
about the origin of the Universe. In particular one expects that the cosmological 
singularity should be avoided. These facts make investigation of the quantum 
gravitational models of the Universe so important and fascinating. 

Many efforts have been made to tackle this issue. One of the most promising 
propositions was the application of the Superstring Theory (ST). In this approach, 
called the pre-Big Bang cosmology [1], @, f|, |4j], the t-duality has been applied 
leading to the avoidance of the initial singularity. The resulting dynamics possess 
a non-singular type of evolution with the minimal universe scales corresponding to 
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the self-dual point of the t-duality transformation. Another interesting approach 
is based on Loop Quantum Gravity (LQG) || which is a promising candidate 
for the quantum theory of gravitation. The methods of LQG applied in the 
cosmological context lead to the avoidance of the initial singularity ||. In this 
scenario the universe is initially in the low energy contracting phase, reaches 
a minimal size and due to quantum effects evolves toward the expanding low 
energy phase |7j ||, |J. In the case with the cosmological constant this model 
was studied in |I(J. It was shown that bouncing behaviour occurs for a dense 
subspace of the physical Hilbert space and semi-classicality is preserved across 
the bounce |TT], The similar quantum bouncing behaviour is also recovered 



in Bojowald's approach, for a recent review of this issue see [I_3|, [14| . 

In this paper we combine these two approaches to construct a non singular 
model of the universe without the need for the non-deterministic t-duality trans- 
formation. This idea, that the LQG effects can have important influence for the 
pre-Big Bang cosmology, was originally proposed by De Risi et al. [15| 

In the present paper we investigate the model of De Risi et al. in more detailed 
way and show some new resulting consequences. Especially we find the analytic 
solution of this model and study numerically effects of the dilaton potential. 

The presented approach is based on replacement of the starting point action. 
Namely the LQG is the background independent quantisation of the Hilbert- 
Einstein action 

= 77- / d 4 x^R. (1.1) 



2k _ 

On the other side, ST of the closed strings leads to the gravitational action in 
the form |, [17] 

1 



R+(V<f ) y-V( ( j ) )--H 2 



;i-2) 



where H 2 = H^ ua H flva and H^ va = 3\d[^B UQ -\. Here beside the metric field 
two additional fields occur. Namely the dilaton scalar field and the asymmetric 
Kalb-Ramond field B^ u . In case when these fields vanish action (|1.2|) take a form 



of the action (|1 . 1| ) . Otherwise the scalar field is non- minimally coupled with the 
curvature and can lead to the important modifications of dynamics. It is worth 
to note that action ( |1.2|) is related with the Brans-Dicke theory with the B-D 
parameter u — — 1. This means that it can be considered as a theory with the 
dynamical coupling constant defined as 

9 2 s = et (1.3) 

1 Early works in LQC also indicate singularity resolution for brane collision model fllq]. 
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In the case when <p — > — oo then g 2 s —>■ and asymptotic freedom occurs. On the 
other side <p —>■ oo leads to the confinement. In the present Universe we expect 
g 2 s ~ const. 

In the presented approach we consider action ( |1.2j) instead of (|1 . 1|) as a start- 
ing point for the loop quantisation. This means that we can consider ST as a 
fundamental theory and the LQG as a background independent method of the 
quantisation of the metric field. However this interpretation is not unique. In 
fact we do not even need to assume ST to be valid but rather to postulate form 
of the action. In this paper we restrict our preliminary investigations to the 
graviton-dilaton sector in the D — 4. However the Kalb-Ramond field is not 
coupled non-minimally to the gravitational field and has less important influence 
for dynamics. 

The organisation of the text is the following. In section |2| we introduce 
holonomy corrections to the action ( |1.2| ) in the Einstein frame. Then in section |3] 
we investigate analytically resulting cosmological model. In section |] we improve 
the analytical model by introduction of the additional potential for the dilaton 
field. In section [| we summarise the results. 



2. String cosmology with holonomy corrections 



In this section we show how to introduce phenomenological effects of the quan- 
tum holonomies in the classical equations of motions. We also give the heuristic 
explanation of the form of these corrections. In result we obtain the phenomeno- 
logical Hamiltonian which contains information about the effective form of the 
quantum gravitational effects. The resulting dynamics trace the mean value of 
the quantum state. However to obtain information about dispersions and approve 
semi-classical considerations full quantum treatment should be performed. 

To introduce loop modifications to the action (|L2|) we choose the strategy 



already applied to the theory with non-minimally coupled scalar fields [IS, 19 



The presented approach is also equivalent to that presented in the paper [|15 
Namely, we perform the conformal transformation in the form 

20 



D-2 



= exp 



which transforms the effective string action (|L 
the Einstein frame where 



(2.1) 



from so called string frame to 



S 



1 



2k d 
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D 
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(2.2) 
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We restrict now our preliminary consideration to the graviton-dilaton sector in 
the D = 4 case, then 



S 



dtL 



2^ / d ^~~ 9 



i?_I ( V0) 2 -r(0)e^ 



(2.3) 



In this case K±g 2 s — k — 8itG where G is the classical Newton coupling constant. 
To quantise this theory in the background independent way we perform the Leg- 
endre transformation of the L to the canonical formulation. Then we express 
the obtained Hamiltonian in terms of the Ashtekar variables (A, E) which 
take value in su(2) and su(2)* algebras respectively. The full Hamiltonian is the 
sum of the gravitational part Tic and dilaton part Ti^. In terms of the Ashtekar 
variables the Hamiltonian for general relativity is a sum of constraints 



He 



d 3 x(N*G t + N a C a + Nh sc ), 



(2.4) 



where 



C a 
G, 



E^ ab 

DaE? 



(2.5) 



and the scalar constraint has a form 



d 3 *N(x)h s 



2^ / d ^ 



y/\detE 



■^ k F k ab - 2(1 + 7 2 



EfE\ 



yf\ det£ 



:KIK 



I-'] 



(2.6) 



with the curvature of the Ashtekar connection F = dA + | [A, A] and where 7 is 
the Barbero-Immirzi parameter. The dilaton Hamiltonian has the form 



d 3 *N(x) - 



where we have defined 




+ Vldet E\U$) ) (2.7) 

(2.8) 
(2.9) 



The next step is to define a spacetime symmetry. It is worth to note that 
if we assume a given symmetry in the string frame it will be the same in the 



- 4 - 



Einstein frame. The only difference is the transformation of time t — ► i and the 
scale factor a — > a for which the form of the metric holds. We choose the flat 
FRW k = spacetime for which metric can be written as 

ds 2 = -N 2 (x)dt 2 + q ab dx a dx b (2.10) 

where N(x) is the lapse function and the spatial part of the metric is expressed 
as 

q ab = Sijuiul = ~a 2 {t)°q ab = a 2 (t)^<4 (2-11) 

In this expression °q ab is the fiducial metric and °u l a are the co-triads dual to the 
triads °ef, °tu l (°ej) = <5] where °uj 1 = °uj l a dx a and °6i = °e°:d a . In this case the 
Ashtekar variables take the form 

A = Y + 1 K = cV ~ 1/3 °uj i a T i dx a , (2.12) 
E = Vldetgle = pV^ 3 ^ (2.13) 

where Vq is the volume of fiducial cell. The volume Vq is just a scaling factor 
and can be chosen arbitrary in the domain Vq £ R + . The physical results do 
not depend on the choice of Vq. The pair (c,p) are canonical variables for the 
gravitational field and can be expressed in terms of the standard FRW variables 
(c, \p\) = ('ydV^ 3 , d 2 V 2 ^ 3 ). Now it is straightforward to calculate the full classical 
Hamiltonian in the canonical variables (c, p, <f), 7r) 

2 

K = -^VW\c 2 + ] ^ 2 + \p\ 3/2 U(4>) (2.14) 

where we assumed homogeneity of the field <ft and chosen gauge N(x) = 1. Loop 
quantisation of such a model in case of the free field was performed in the works 
0, p|, P . Here we only sketch the main steps and for the detailed considerations 
we send the reader to the mentioned papers. To quantise this theory in the 
background independent way one introduces holonomies of connection A 

h a [A}=Vexp / A where 1-form A = TiA i a dx a (2.15) 

J a 

and conjugated fluxes 

F l s [E]= [ dF i where 2-form dF { = e abc E?dx b A dx c (2.16) 
Js 
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where a.SGE and lir^ = a,- where a,- are the Pauli matrices. From this definition 



we can calculate holonomy in the direction °e\d a and the length /j,V 



1/3 



h 



(p) 



exp 



1/3 



r iC y - 1/3 Xdx a = I cos (^) + 2r, sin 



(2.17) 



where we used the definition of the Ashtekar variable A ( |2.12| ). From such a 
particular holonomies we can construct holonomy along the closed curve a = djj . 
This holonomy can be written as 



1 — Hj L J L J 



(2.18) 



Now it is straightforward to show that the field strength can be expressed as 



tr 



ab 



-2 lim ■ 

Ar-fO 



2/3 



J V V 



b- 



The trace in this equation can be explicitly calculated 



tr 



□i,- 



I 



^kij 



sm 



(/ic) . 



;2.i9) 



(2.20) 



In Loop Quantum Cosmology the limit Ar — > in the formula (|2.19|) does not 
exist because of existence of the area gap. The area gap corresponds to the 
minimal quanta of area A = 2^/3^1^ PH] . So instead of the limit in equation 
( |2.19| ) we should stop shrinking the loop at the appropriate minimal area A. 
This area corresponds to the area Hj intersected by the loop and to take account 
discreetness of the space we should perform a limit EL- — > A. Now, we must 
connect the area Elj with the length /i of the loop edge. We can choose that area 



correspond to the physical area a 2 /i 2 . 
we have in the limit 



So in this case EL = V r 



2 /3~2, ,2 







a yU 



v = v '■= \ n 




(2.21) 



We will use this value in the further consideration. However we must to mention 
that the choice of the function \x is not unique and leads to ambiguities. However 
function /2 used here leads to the proper classical limit and follows from loop 
quantization 0. Special situation corresponding to so called \i — scheme where 
\i = Ho =const. In such a case /i depend explicitly on the factor Vq. However 
it was shown that fi — scheme does not have the classical limit and cannot be 
treated physically P, In general ambiguities of fi(p) come from considera- 



tion of the so called lattice models [23| for which in general p(p) oc \p\ n where 
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n G [—1/2,0] Ul3j . However recent results indicate that the only p, approach is 
consistent quantization scheme [ 24j] . The issue of other kind of ambiguities was 
addressed in the papers p5| ^6 



Based on the above considerations we conclude that phenomenological effects 
of the quantum holonomies can be introduced by a replacement 

sin (p,c) 



2.22 



in the classical expression. In fact another type of corrections occur in the Loop 
Quantum Cosmology. However there is no corrections in the flat models. Namely 
the inverse volume corrections which should be in general included have shown 
to be physically inconsistent in the context of the flat models. This is the result 
of the non-compact and flat manifold S. In the cosmological models with the 
internal curvature like FRW k = ±1 inverse volume corrections should be also 
taken into account |27 |. 



3. Non-singular universe from vacuum 

The introduction of the holonomy corrections investigated in the previous section 
lead to the effective Hamiltonian in the form 



n 



eff 



K 4 T 



sin (/xc) 



2|p| 3/2 



(3.1) 



where we assumed homogeneity of the dilaton field. This Hamiltonian fulfils 
constraint 7Y c ff = 0. The equations of motion can be now derived with use of the 
Hamiltonian equation / = {/, 7i c s} and for the canonical variables are 



P = {P, ^eflf} 
C = {C, Hes} 



1 V 



sin (fj,c) cos (jic) 



sin {lie) 



- sgn(p) 



k 4 7 



7T- 



IpI 



5/2 • 



eff 



\P\ 



-3/2 



71"; 



7T; 



(3.2) 

(3.3) 

(3.4) 
(3.5) 



where all derivatives are in respect to the time t. The exact solution for the 

1/3 



variable p can be find [3^] and has a form 
pit) = sgn(p) 



(3.6) 
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This corresponds to the bouncing type behaviour in the Einstein frame. In Fig. 
|l] we show evolution in this frame for the scale factor a and the Hubble parameter 
H 



1 da 
a dt ' 




time in the E-frame 



Figure 1: The top curve (red) presents the bouncing solution for the a(i). Bottom 
curve (blue) corresponds to the time i evolution of the Hubble parameter H = A % . 



The scale factors and times in the Einstein and string frames are related by 

a = e 0/2 a, (3.7) 
dt = e H2 dt. (3.8) 



Therefore to analyse the evolution in the string frame the solution for the dilaton 
field is required. This can be obtained analytically by the simple integration 



Til 



M 3/2 

dt 



i 7 2 A + 1 2 ^ 



sgn(7T7)arcsh 



3t 



(3.9) 



The solution depends on sgn(7r^) leading to monotonically increasing or decreas- 
ing functions. In the case tt^ < there is no birth of the classical universe. 
Therefore in the further considerations we choose tt^ > 0. In this case <p(t) is a 
monotonically increasing function. We show this behaviour together with evolu- 
tion of the energy density 

12 



(3.10) 
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in Fig. |2|. The important feature of this evolution is that energy density does not 
diverge. The maximal value of the energy density 

<°c = — W (3-11) 

corresponds to the hot Big Bang which is here the transitional non-singular phase. 
This result is with agreement with the results of the paper [TTH]. 




time in the E-frame 



Figure 2: Monotonic (red) curve presents the evolution of the dilaton field 4>{t). Picked 
(blue) curve shows the time i dependence of the dilaton energy density p = (ft 2 /2. 

The subsequent quantity that we can investigate is time t dependence of the 



coupling constant ( |1.3|) which is given by 

^ = eXP {^ Sgn( ^ )arCSh (7^) 



(3.12) 



This function rapidly decreases to zero for negative value of time t leading to the 
asymptotic lack of the interactions in the pre-Big Bang phase. In the post Big 
Bang phase the coupling constant g 2 grows monotonically. This is rather non 
required behaviour, because we rather expect the coupling constant to saturate. 
In the next section we show that when the dilaton potential is present such 
a scenario can occur. The presented here model can be rather treated as an 
analytical toy model which requires some further modifications and analysis. 
The evolution of the scale factor in the string frame takes a form 



V Q exp 



sgn(7Tj 



-arcsh 



3f 



7 



.9 a " 9 "*? 



1/6 



(3.13) 
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We show this behaviour together with g 2 (t) in Fig. In the early times t — > — oo 
the scale factor a(t) goes to zero. This state can be interpreted as a vacuum state 
without interaction and with vanishing energy density. It is important to note 
that this is the classical state and the quantum effects are important only in the 
neighbourhood of the Big Bang transition state. 




time in the E-frame 



Figure 3: The top curve (blue) shows the evolution of the scale factor a(t). The 
bottom curve (red) shows the time t dependence of the coupling constant g 2 . 



All above considerations we have performed in the time t. However we can 
back to the time t by integration 



t 



V3 



sgnlrr , 



sgn(7r^ 
exp < arcsh 



V3 



3t 



7VA 



(3.14) 



This is the monotonic function of t so the results obtained with time t are equiva- 
lent with these in time t. Considering asymptotic behaviour of the formula ( ^.13 ) 
and ( |3.14j ) we recover evolution of the post-Big Bang Universe 



,0.577 



a(t) ocfVsRif 

for t ^> 0. The asymptotic behaviour of the dilaton has a form 

(j){t) oc In (t) . 



(3.15) 



(3.16) 
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In the paper |l5j it was suggested that thanks to LQG effects, the dilaton in 



the string frame can approach a constant value for late times. However solution 
( |3.16| ) indicates that such a scenario is not realised in this simple model. The 
increase of dilaton field for the late times is however logarithmically slow. For 
the late times the coupling constant behave as 

g] oc t^- 1 « t°' 732 . (3.17) 

In the next section we show how to improve this model to obtain the satura- 
tion of the coupling constant for the late times. 



4. Possible evolutionary scenario: Prom asymptotically free 
to dark energy dominated universe 

In the previous section we have considered the model with the free dilaton field. 
Resulting dynamics was described by a simple analytical solution. To generalise 
this model it is naturally to add a non-vanishing potential function. This makes 
dynamics more complicated to analysis and requires numerical methods to be 
used. In this section we restrict ourselves to study qualitative modifications in 
dynamics that can occur when dilaton potential is present. 

The simplest form of the potential energy that we can assume is the mass 
potential 

2 

TY) 

V(<j ) ) = —(<f ) - ( f ) o) 2 + Vo (4.1) 

or the 4 type interacting potential. In fact these two possibilities lead to the 
same qualitative prediction and we can concentrate on the potential ( |4.1| ). In this 
case potential U(4>) takes a form 

2 

U{<t>) = —(0 - 0o) 2 e^ + Voe^*. (4.2) 

This potential together with the possible evolution path of the field is sketched 
in Fig. f|. The dynamics of the field <fi is governed by the equation 

d<j> 

This equation can in general contain quantum corrections form the inverse volume 
operator. However they are physically inconsistent in the case of the flat models 



(4.3) 




and are not expected to give contribution in this peculiar case. The dynamics of 
the dilaton filed is therefore governed by the classical formula. 

The potential U(4>) has two local minima, for cf) = <fto and for <fi — ► — oo. 
Lets imagine scenario in which field starts its evolution with asymptotically free 
state cf) — > — oo, g 2 s — > 0. In this phase the scale factor a — > and energy 
density p — > 0. We can call this state asymptotically free vacuum. In this phase 
there is no interactions between particles and zero energy density in agreement 



with the hypothesis of "asymptotic past triviality" [ 29| . In the same time in 
the Einstein frame the scale factor decreases leading to the negative value of the 
Hubble parameter H. This leads to the negative friction term in equation QL3|) 
resulting motion of the field toward higher values. The dynamics of the scalar field 
can be then approximated by + 3H<p pa 0. Then energy density increases and 
the dilaton field reaches the top of the potential hill and subsequently falls down 
to the potential minima. This stage corresponds to the hot Big Bang phase when 
energy density reaches its maximum. However there is no singular behaviour 
here. When the dilaton field falls to the the potential well it dumps because H 
starts to be positive. It means that the non-singular hot Big Bang corresponds 
to the bounce in the Einstein frame. 

Let us investigate now the case with Vq = and <fto = 0. The local maxima 



-12- 



for such a potential is placed in 



what gives 



Ufa 



2 

m 2 1 
k 4 e 2 



Because for physical evolution p < p c , we obtain relation 

-,2 



m 2 < 



(4.4) 



(4.5) 



(4.6) 



27r 7 3 G , 4 

where K4 = 87^4. This relation must be fulfilled to permit the transition of the 
field <j) to the potential well around 
it gives 



m 2 < 



- 0. In particular case 7m = 0.12738 |28 
^ (4.7) 



0= 


0) 


0max) 


a(t = 


0) 


= 1.0, 


P (i= 


0) 


= Pc 


t(t = 


0) 


= 0.0. 



Now we choose m = 25/y/Gl and boundary conditions in the form 

(4.8) 
(4.9) 
(4.10) 
(4.11) 

With such a setup we perform numerical integration of the equations of motion. In 
Fig. [5| we show calculated evolution of the scale factor and the Hubble parameter 
in the Einstein frame. This behaviour is qualitatively similar to the earlier studied 
case without the potential. However here the symmetry t —t is broken. In the 
left panel in Fig. |^ we show dependence between time in the string frame and time 
in the Einstein frame. We see that for the positive times this dependence is one 
to one. This is the result of the dilaton stabilisation which gives e^/ 2 — > 1 in the 
transformation ( |3.8j ). This fact can be seen from the right panel in Fig. |B| where 
evolution of the dilaton field is shown. Stabilisation of the dilaton field lead also 
to the saturation of the coupling constant g 2 s as it was in the left panel in Fig. |7|. 
In the right panel in Fig. [7] we show corresponding evolution of the scale factor 
in the string frame. In the small time scales this evolution undergoes oscillations 
resulting from the dilaton field damping oscillations. In average scale factor 
behaves as for the universe dominated by dust with the equation of state p = wp 
where (w) = 0. The left panel in Fig. show the evolution of the parameter w. 
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time in the E-frame time in the E-frame 



Figure 5: Left : Evolution of the scale factor in the Einstein frame. Right : Evolution 
of the Hubble parameter in the Einstein frame . 




-3 -2 -1 1 2 3 0.5 1 1.5 2 2.5 3 

time in the E-frame 1 

Figure 6: Left : Dependence between time in the S-frame and time in the E-frame. 
Right : Evolution of the field <p in the time t. 
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Figure 7: Left : Temporal evolution of the coupling constant g 2 a . Right : Evolution 
of the scale factor in the string frame. 



- 14 - 




Figure 8: Left : Evolution of the w parameter in the equation of state p = wp. Right 
: Evolution of the energy density of the dilaton field. 

In the right panel in Fig. |8| we show evolution of the dilaton energy density p. In 
the model considered above we have obtained expected saturation of the coupling 
constant g% — > const and classical post-Big Bang state with (w) = 0. Now we will 
add non-vanishing constant contribution Vq 7^ to the dilaton potential. Such 
a modification is important only close to the local minimum point = what 
correspond to the post-Big-Band phase. In the bottom of the potential well the 
dilaton field behaves like dark energy with the equation of state p = wp where 



w 



L. - utf) -Ufa) 



-Wo) 



Vn 



(4.12) 



leading to the classical de Sitter phase. This scenario we study numerically in 
the case Vo = 0.01/Crf. In the left panel in Fig. |] we show resulting behaviour 
of the parameter to, which is damped to the constant value w — — 1. In the right 
panel in Fig. |9| we show corresponding evolution of the scale factor in the string 
frame. The dilaton field can therefore behave like inflaton field or dark energy. 



5. Summary 

In this paper we have analysed the low energy superstring theory with Loop 
Quantum Gravity effects as the tool to investigate the early Universe. This 
approach was initially proposed by de Rissi et al. |T5| . Based on this approach 



we had derived the cosmological model and solved it analytically. The important 
result which occurs is that an initial singularity state is avoided without need for 
the t-duality transformation. In this scenario the universe starts its evolution from 
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Figure 9: Left : Evolution of the w parameter in the equation of state p = wp in the 
case Vo = 0.01/G 2 . Right : Corresponding evolution of the scale factor in the string 
frame. 

asymptotically free vacuum. Namely, the state characterised by the vanishing 
coupling constant g 2 s — > 0, energy density p — » and the scale factor a — > 
for t — > — oo. This picture is the realisation of the "asymptotic past triviality" 
hypothesis. The universe starts its evolution from such a state and then evolves 
toward the non-singular Big Bang phase and subsequently to the current state. 
However in this simple model the coupling constant monotonically grows during 
evolution. The effective evolution is asymptotically described by a(t) oc t 1 ^ for 
t — > oo. 

We have also proposed improvement of this model introducing potential for 
the dilaton field. In this case we also start from the initial asymptotically free 
vacuum. However in the further stages of the evolution, the dilaton field is 
constrained by the potential part e^. After the non-singular Big Bang phase, 
characterised by the maximal and bounded energy density, the dilaton field ends 
its evolution in the bottom of potential well. This leads to the saturation of the 
coupling constant gl — *> const. The dilaton in this phase behaves like dark energy 
with the equation of state p = —p. This leads to the emergence of the de Sitter 
phase which can explain present cosmological observations. 
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